Parity Nonconservation in Odd-isotopes of Single Trapped Atomic Ions 
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We have estimated the size of the light-shifts due to parity nonconservation (PNC) interactions 
in different isotopes of Ba + and Ra + ions based on the work of Fortson [Phys. Rev. Lett. 70, 
2383 (1993)]. We have used the nuclear spin independent (NSI) amplitudes calculated earlier by us 
[Phys. Rev. Lett. 96, 163003 (2006); Phys. Rev. A 78, 050501(R) (2008)] and we have employed 
the third order many-body perturbation theory (MBPT(3)) in this work to estimate the nuclear 
spin dependent (NSD) amplitudes in these ions. Ra + is found to be more favourable than Ba + for 
measuring both the NSI and NSD PNC observables. 

PACS numbers: 



Parity nonconservation (PNC) in an atom arises 
mainly due to the exchange of the Z° boson between 
the electrons and the nucleus The interaction lead- 
ing to such an effect consists of two parts. One of them 
is nuclear spin-independent (NSI) and the other is nu- 
clear spin-dependent (NSD). In addition, the interaction 
between the electrons and the nuclear anapole moment 
(NAM) in an atomic system which is NSD in charac- 
ter, can also give rise to PNC [HHf. Wood et. al. 0] 
have reported the observation of the NAM in atomic Cs. 
For an atomic system, the contribution of the NSI and 
NSD components can be extracted by using the addition 
and subtraction of two separate measurements respec- 
tively. However, it is not possible to distinguish the NSD 
contributions from the NAM or Z° exchange from these 
measurements, even though the latter is typically larger 
than the former in heavy atoms. Henceforth in this work, 
the term NSD interaction will refer to the combination 
of these two interactions. 

Although Wood et. al. [2| claim to have observed 
the NAM in atomic Cs, it has so far not been possible 
to explain this observation from well established nuclear 
data 3. Therefore, there is an urgent need to study the 
NSD part of atomic PNC. In atomic transitions when 
both the NSI and NSD interactons contribute simultane- 
ously, the typical strength of the NSD component is much 
smaller than its NSI counterpart. The overall magnitudes 
of the both interactions grow rapidly with the atomic size 
[H . It may therefore appear that heavier atomic systems 
would be a natural choice for studying the PNC transi- 



tions. But in a certain cases, the PNC observables can 
also be enhanced due to degenracies To date all 

the observations of atomic PNC have been reported in 
neutral systems However, Fortson in 1993 proposed 
that a single trapped and laser cooled ion can be used 
for measuring PNC with an accuracy that is comparable 
to that of their neutral counterparts 5]. This proposal 
is based on the interference of the PNC induced electric 
dipole amplitude {Elp^c) with the electric quadrupole 
(E2) amplitude in the 6s 2 Si/ 2 —> 5d 2 D 3 / 2 transition in 
Ba + which leads to a PNC induced light shift. This ap- 
proach is being pursued at KVI [fj, 0] in an experiment 
with the ion of the next heavier alkaline earth element 
Ra + , which has low- lying optical transitions. This ef- 
fort has been preceded by many recent experimental and 
theoretical studies on properties related to PNC (6- 12|. 
and it appears that a high precision result for PNC in 
these ions is possible. In addition, a novel idea of observ- 



ing only the NSD PNC interaction in the S1/2 —> D. 
transitions of these ions has been reported in Ref. [1 
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In this paper, we compare the possible light-shifts due 
to PNC in the 2 S 1/2 ^ 2 D 3/2 and 2 S 1/2 -> 2 D 5/2 tran- 
sitions for Ba + and Ra + using our calculated results re- 
ported earlier [^, IB - 12] as well as from this work. Only 
the NSD PNC amplitudes are evaluated here using the 
third order many-body perturbation theory (MBPT(3)). 
The results of these studies will aid in the choice of the 
hyperfine level transitions involving the S1/2 and -03/2,5/2 
states in different isotopes for performing PNC experi- 
ments in order to extract the NSD PNC observables in a 
systematic and unambiguous way in the two above men- 
tioned ions. 
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To measure the NSD PNC contribution, we propose to 
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carry out measurements on 135,137 Ba (/ = 3/2), 139 Ba 
(I = 7/2), 225 Ra (/ = 1/2), 223 Ra (/ = 3/2) and 229 Ra 
(/ — 5/2). All these isotopes have long nuclear lifetimes 
except 229 Ra, where the measurements can be performed 
in an on-line facility. 

In a single ion experiment, the observable is the PNC 
electric dipole transition induced ac-Stark shift i.e. the 
PNC light shift in the ground state Zeeman sub-levels 
when a laser drives a E2 transition between the ground 
state and a meta-stable state. Generalizing the expres- 
sion given by Fortson Q, the measured PNC light shift 
in the S\/ 2 — > ^3/2,5/2 transitions in a single ion experi- 
ment will have a small shift due to PNC as 
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and a much larger shift due to the quadrupole coupling 
between the same two levels, 
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for a given sub-level mj? of a hyperfine state of angular 
momentum F, where lvq and ui are the atomic and opti- 
cal frequencies, respectively. The Rabi frequency for the 
PNC-induced-dipole transition is given by 
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and the Rabi frequency for the quadrupole transition is 
given by 
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where £ is the applied electric field. We use atomic unit 
(au) through out this work unless mentioned explicitly. 

To the best of our knowledge, the E1 PNC ampli- 
tudes for the 5*1/2 — > D3/2 transitions due to the NSD 
interaction for these ions have not been calculated so 
far. Preliminary calculations of these amplitudes for the 
S1/2 — > D5/2 transitions are given in Ref. [3] using the 
rclativistic configuration interaction (CI) method, but 
the expression used in that work for the NSD-interaction 
is not compatible with our analysis. In what follows, we 
employ a relativistic MBPT(3) method (as described be- 
low) here to evaluate these amplitudes in a systematic 
fashion for both the NSD PNC transition amplitudes. 

The NSD part of the PNC interaction Hamiltonian is 
given by 
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where Gf is the Fermi constant, p nuc {r) is the nuclear 
potential, a is the Dirac matrix, I is the nuclear spin and 



lZ a is a dimensionless constant which has information 
about NAM. The Elp^c amplitude due to this Hamil- 
tonian can be written as 
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where the subscript and 1 represent the atomic wave- 
functions and their first order corrections due to Hpj^®, 
i and / represent the valence orbitals in the initial and 
final states, respectively, and D is the electric dipole op- 
erator. 

In the MBPT method, we define wave operators il v .o 
and fi^i to calculate the unperturbed fl^w )) and per- 
turbed )) wavefunctions as 

00 
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where is the Dirac-Fock (DF) wavefunction obtained 
using the Dirac-Coulomb Hamiltonian. 

We use the generalized Bloch equation given below to 
calculate the unperturbed wavefunctions in our MBPT 
formulation (for n > 1) 



m— 1 

where Hq is the DF Hamiltonian, V r is the residual 
Coulomb interaction, f2^°g = 1 and P and Q are the 
projection operators in the model and orthogonal spaces, 
respectively; i.e. 

P =|$„)($„ I and Q=l-P. (9) 

Following the similar procedure (for n > 1), we get 
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where = QHgfigP. For the MBPT(3) approxima- 
tion, we consider terms up to n = 2 to evaluate -Elpjvc 
given by 

ElpNC = (i-i-) 
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In the angular momentum relations, we express 

EIpnc = ((JuI^FiMilDlff + D 2 eff \(J f ,I);F f Mf) 
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where F is the total angular momentum due to the elec- 
tron angular momentum (</) and the nuclear spin (I) 
with its azimuthal component M. Here q ~ —1,0,1 



J 



depending upon the values of Mi and Mf. M. = 
((Ji,7);F i ||£»i // + £»2 // ||(J /) /) ; J? / ) is the reduced ma- 
trix element of effective rank one operators with 
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and D 2 eff = ^ ™. (13) 



The above expression is non-zero for Fi = Ff , Ft ± 1 . 
With the expansion H^g = it gives 



and 
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where r\ — \JI(I + 1)(2J + l)(2Fi + 1)(2F/ + \)/norm, ej represents orbital energy for i and the matrix element in 
terms of single particle orbitals is given by 
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for P(r) and Q(r) being the large and small radial com- 
ponents of Dirac wavefunction and a is the Pauli spinor 
with (m\ \a\ I /c/) = (ftj + ft/ — l)(Kj||C 1 | |k/) for the Racah 
tensor C and the relativistic quantum number ft. 

In the method proposed by Fortson, the PNC electric 
dipole transition induced light shift is observed as a fre- 
quency shift of the ground state Larmor frequency in spin 
zero isotopes. Any fluctuation of the quadrupole transi- 
tion induced light shift does not modify the uncertainty 
of the PNC light shift measurement, since it cancels out 
in the Zeeman transition due to its tuf dependence (see 
Eqs. {1} and ([2])). However, in the case of non-zero spin 
isotopes of these ions, the ±rriF transitions are forbidden 
by the E2 selection rules. 

The estimated magnitudes of the PNC light shifts for 
both NSI and NSD in the Si/ 2 — > 03/2,5/2 transitions are 
given in Table U (only non negligible values are given). 
For the calculations, we have used our earlier results for 
NSI PNC and electric quadrupole transition amplitudes 
and the NSD amplitudes have been evaluated 
in this work. On practical grounds, we have taken the 
strength of the electric field as 2 x 10 6 V/m; this value 
optimizes the quenching rate for Ba + , and R a ~ 0.2, 
obtained from the Cs NSD PNC studies 

It is clear from the above table, the NSD PNC light 
shift in the case of the S1/2 —> D3/2 transition is about 
a few mHz (see Table [T| which requires the measure- 
ment of the combined PNC light shift to less than 0.5% 



precision. The allowed spin flip transitions in these iso- 
topes are associated with not only the desired differen- 
tial PNC light shift but also a much larger differential 
quadrupole light shift. Thus to achieve a precision be- 
low 0.5%, all the sources of temporal variation of the 
quadrupole light shift need to be stable with uncertainty 
well below a percent. Only in such cases the measure- 
ments of the quadrupole and PNC light shifts within a 
short time interval can provide NSD PNC light shifts 
with the desired accuracies. However, a suitable choice 
of the hyperfine states and the Zeeman sublevels allows 
to avoid the systematics from the quadrupole light shift. 
In the F = 2(Si/ 2 ) — > F' = 3(D 3 / 2 ) transition in spin 
I = 3/2 isotopes, the sublevels tuf = 1, of F = 2(Si/ 2 ) 
will have the same quadrupole light shift and this will 
reduce the systematic error in the NSD PNC measure- 
ment. An unambiguous measurement of NSD PNC will 
indeed be a challenge, since it is necessary to know the 
NSI PNC part in the same isotope. An alternative is to 
consider the F — 3(^/2) — > F' = 2(D 3 / 2 ) transitions in 
139 Ba+ (I = 7/2) and 229 Ra+ (I = 5/2). In these iso- 
topes rap = 3,1 sublevels (F — 3,S 1 / 2 ) experience the 
same quadrupole light shift while mp = 2 is free from 
the quadrupole and PNC light shifts. Thus it would be 
possible to extract both contributions to the PNC by 
driving the spin flip transitions after preparing the ion in 
the tuf = 2 state. 

In principle an unambiguous measurement of NSD 
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TABLE I: Induced light-shifts due to E2 and PNC interactions in different isotopes of Ba + and Ra + . 1Z a ~ 0.2 from the Cs 
data 2, 3] and electric fields as 2 x I0 6 V/m are considered for the practical estimations. 

Transition lj~l/2\El^\j^l/2) {ji\\E2\\j f ) F, F f {Fi\\El^\\F f ) m F Aoj^ uad /2-K Acj££ j c /2tt Aw^/2tt 
ji->jf (xlQ-^jeap) [6, .10] (au) [11, .12] (x lQ- 14 ieao/^ a ) (MHz) (Hz) xlQ- 4 (Hz) 



'/ 137 Ba+ (/ = 3/2) 



6s — y 5dq/o 




2.46 


12.74 
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3 


97.18 


1 


7.62 


-0.24 


6.08 








1 


2 


-89.40 


1 


9.53 


-0.37 


-14.03 


6s — > 5d 5 / 2 







15.96 


2 


3 


1.37 


1 


11.13 





0.09 


139 Ba+ (/ = 


7/2) 




















6s — > 5dQ/o 




2.46 


12.74 


3 


3 


103.57 


1 


5.07 


0.42 


12.29 
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2 


-105.55 


3 


3.6 


-0.47 


-14.75 
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3.6 


-0.1 


2.95 


6s -> 5d 5/2 
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0.608 
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12.87 





0.08 
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12.87 





-0.02 


225 Ra+ (7 = 


1/2) 




















7s -> 6d 3/2 




46.4 


-14.87 


1 


2 


991.75 


1 


38.95 


9.97 


155.6 


223 Ra+ (7 = 


3/2) 




















7s ->■ 6d 3/2 




46.4 


-14.87 


2 


3 


1173.45 


1 


22.03 


-4.7 


73.35 








1 


2 


-1017.49 


1 


27.54 


-7.05 


-159.64 


7s ->■ 6d 5/2 







-19.04 


2 


3 


-17.52 


1 


32.08 





-1.1 


229 Ra+ (7 = 


5/2) 




















7s ->• 6d 3/2 




46.4 


-14.87 


3 


2 


325.94 


3 


29.44 


-4.19 


45.56 














1 


29.44 


0.83 


-9.11 










2 


2 


1147.73 


1 


17.42 


7.42 


152.18 


7s ->■ 6d 5/2 







-19.04 


3 


2 


-5.12 


3 


17.5 





-0.72 














1 


17.5 





0.14 



PNC would be possible in the Si/ 2 — > -D5/2 transi- 
tion, but the size of the observable light shift is very 
small; hence its realization is a challange. In Table U 
the light shifts in the S1/2 — > D 5 / 2 transitions with dif- 
ferent hyperfine states are given for various isotopes of 
these ions where it is feasible to drive spin flip transition 
between the magnetic sublevels avoiding the differential 
quadrupole light shift. It shows that an experimental 
uncertainty below 0.05% of the Larmor frequency shift 
measurement is essential to extract NAM result. 

In conclusion, we have given the estimated values of 
the light shifts due to the nuclear spin independent and 
dependent parity nonconserving interactions for various 
isotopes of singly ionized barium and radium isotopes. 
We have shown that the size of these effects would be 
rather small for Ba + , and therefore it would be quite 
challenging to observe the nuclear spin dependent par- 
ity nonconserving effect in this ion. But in the case of 
Ra + , the prospects for observing the parity nonconserv- 
ing light-shift due to the nuclear spin dependent inter- 
action are much better. An unambiguous observation of 
light-shift due to nuclear spin-dependent interaction in 
the 7s 2 5*1/2 6d 2 D 5 / 2 transition might be feasible in 
singly ionized radium. Our analysis highlights the iso- 
topes of Ra + and the transitions in them that could be 
suitable for observing the nuclear spin-dependent parity 
nonconservation. 
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